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We systematically reconsider, within an improved "analytic if-matrix model", the extraction of the a = /o(600) and /o(980) 
masses, widths and hadronic couplings using new Kci = K — ► -K%ev e data on 7T7t phase shift below 390 MeV and different sets 
of 7T7T — > tttt/KK scatterings data from 400 MeV to 1.4 GeV. Our results are summarized in Tables [1] [2] and [5] In units of 
MeV, the complex poles are: M a = 452(12) — i 260(15) and Mj = 981(34) — i 18(11), which are comparable with some recent 
high-precision determinations and with PDG values. Besides some other results, we find: \g a K+K~ \/\9<jit+tt~ I = 0.37(6) which 
confirms a sizeable g a K+K- coupling found earlier, and which disfavours a large 7T7t molecule or four-quark component of the a, 
while its broad 7T7t width (relative to the one of the p-meson) cannot be explained within a qq scenario. The narrow 7T7t width of 
the /o(980) and the large value: \gtjc+K- l/lff/ir+n— I ~ 2.59(1.34), excludes its pure (mm + dd) content. A significant gluonium 
component eventually mixed with qq appears to be necessary for evading the previous difficulties. 



1. Introduction 

Understanding the nature of scalar mesons in terms of 
quark and gluon constituents is a long standing puz- 
zle in QCD |1I2) . The problem here is that some 
states are very broad (cr and k [3] (if confirmed) mesons) 
and others are close to an inelastic threshold (/o(980), 
a (980)), which makes their interpretation difficult. Be- 
sides the interpretation within a qq model [1I4I5I6I7I8I9] 
or unitarized quark model [10111) . also the possibility 
of tetraquark states |12ll3ll4ll5ll6j (and some other re- 
lated scenarios: meson- meson molecules |17|18) . meson 
exchange [19]) is considered. In addition, a gluonic me- 
son is expected back to 1975 [20] as a consequence of 
the QCD confinement. An indication of the existence 
of scalar glueballs [212 1] comes from phenomenologi- 
cal studies |22I23I24[ , lattice QCD [25I26J . strong cou- 
pling approach in a Nambu-Jona-Lasinio like-model [27) . 
QCD spectral sum rules (QSSR) |6l9l28l29l30l3"TI52l ?] 
a la SVZ [33134] . some low-energy theorems (LET) 
[32135136] . AdS/QCD [37] and large N c [55159] . Such a 
state could mix with the other qq mesons [5I6I8I40I41] , 
Among the light particles, the a = /o(600) (hereafter 
called cr) meson could be such a gluonic resonance. The 
cr can manifest itself in some effective linear sigma mod- 
els [42143144] or contribute to the low-energy constants 
at 0(p 4 ) of the QCD effective chiral Lagrangian [45], 
while its role in nuclear matter (e.g. nuclear potential) 
from its coupling to nucleons is essential [46] . 
One might expect that the hadronic and two photons 
couplings of these mesons could provide an important 
information about their intrinsic composite structure. 
Indeed, recent analyses of 77 — > tttt data [47] and of 
tttt — > tttt/KK [48] indicate that the a meson could be 
such a gluonic resonance (gluonium/glueball). 
In this letter, we pursue the test of the nature of the 
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cr and of the /o(980) by reconsidering the extraction 
of their parameters (masses and couplings) from tt + tt~ 
scatterings using recent precise data from NA48/2 on 
Ke4 = K — > nireve for the 7T7r-phase shift below 390 
MeV [49] and different sets of 7T7t — > tttt/KK data above 
400 MeV [50151152153154] . 

2. The analytic K-matrix model for 7T7t — > tttt/KK 
In this approach, the strong processes are described by 
a K matrix model representing the amplitudes by a set 
of resonance poles [55]0 . In that case, the dispersion 
relations in the multi-channel case can be solved explic- 
itly, which is not possible otherwise. The model can 
be reproduced by a set of Feynman diagrams, including 
resonance (bare) couplings to 7T7t and KK and (in the 
original model [55 ) 4-point 7T7t and KK interaction ver- 
tices which we shall omit for simplicity in [47] and here. 
A subclass of bubble pion loop diagrams including res- 
onance poles in the s-channel are resummed (unitarized 
Born). In this letter, we discuss the approach for the 
case of : 1 channel 1 "bare" resonance (K-matrix 
pole) and 2 channels © 2 "bare" resonances and we re- 
strict to the SU(3) symmetric shape function. 
In the present analysis, the introduction of a real ana- 
lytic form factor shape function, which takes explicitly 
into account left-handed cut singularities for the strong 
interaction amplitude, allows a more flexible parametri- 
sation of the 7T7t —¥ tttt/KK data. In our low energy 
approach, it can be conveniently approximated by: 

fp(s) = , P = tt,K, (1) 

S + <JDP 

which multiplies the scalar meson couplings to tttt / KK . 
In this form, the shape function allows for an Adler zero 
at s — sap and a pole at odp > simulating the left 
hand cut. 

1 channel © 1 "bare" resonance 

4 Some applications of the model have been discussed in [56 57 . 
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Let's first illustrate the method in this simple case. The 
unitary PP amplitude is then written as: 

G P fp(s) _ G P f P ( S ) 



and the partial propagators 



Tpp(s) 



sr - s 

where Tp p — 
(l-Amp/s) 1 / 2 ; 



(2) 



-G P fp(s) V P (s) ' 

e lbp sin Sp/pp(s) with pp(s) = 
Gp = 9 Per b are the bare coupling 
squared and : 

Im V P = Im (-Gpfp) = -{9pp)G P f P , (3) 

with: (9pp)(s) = below and (6pp)(s) — pp(s) above 
threshold s = Amp. The "physical" couplings are de- 
fined from the residues, with the normalization: 

g 2 Pa = glpp/(l^) . (4) 

The amplitude near the pole so where T>p(so) = and 
V P (s) ^ W p(s )(s - s ) is: 

Gpfp(s a ) 



Tpp(s) 



9p 



9 2 p 



(5) 



so - s -£>'(so) 
The real part of T>p is obtained from a dispersion rela- 
tion with subtraction at s — and one obtains: 
2 



fp(s) 



hp(s) -h P {0) 



(G) 



hp{s) = fp{s)L s i(s)-(a N p/(s+<T D p))L s i(-<7pp), cfnp 
is the residue of fp(s) at —a up and: L sl (s) = 
[ (s — Amp) /mp]Li(s,mp) with Li from [55] . 
2 channels 2 "bare" resonances 
The generalization to this case is conceptually straight- 
forward though cumbersome. Let us consider two 2- 
body channels coupled to 2 "bare" resonances labelled 
a and 6, with bare masses squared sp a and sp b : 

- Let /„(s) , Ub(s) jKa{s) , f Kb{s) be four shape func- 
tions, real analytic in the s-plane, with left cut, and 

Uaa(s), fwbb(s) , fwab(s), fKaa{s), f~Kbb(s) , /if a b(s), 

six functions, real analytic in the s-plane, with right cut. 
Their imaginary parts on the cut for s > Amp are: 

lmU aa (s + ie) = (6p«fl a ){s), (7) 
Iin/U fc (s + it) = (6pn fl h )(s) , 
lmf„ab(s + ie) = (Op^ f„ a / 7rh )(s) , 

and analogous for the 2nd KK channel. 

- Let's define the bare inverse propagators: 

D a (s) = (s Ra - s); D b (s) = (s Rb - s) , (8) 

and the "bare" couplings g^ a , g-Kb, gKa, gxb of the res- 
onances to the channels, through the pure 1-resonancc 
inverse propagators: 

T> a (s) = D a (s) - gl a f^aa{s) ~ g\a f K aa{s) , 

V b (s) = D b (s) - gl b Ubb(s) - g 2 Kb f K bb(s) ■ (9) 

- Let's define the full denominator function T>(s), ana- 
lytic in the s-plane, with right cut s > Am 2 : 



V(s) = V a V b - {gTragrrbfirab " 
= D a P > b — Da {g^b fnbb 



gKagKbfKab) 2 
- 9Kb fKbb) 



Db (Siixa fvaa + 9Ka fKaa) + (fiVfc U 



iraa 
1 

7T a 



1,1, 



+9 Kb fKbb)(g: 



+ 9Ka fKaa) 
9K a gKbfKab) 2 , 



(10) 



p _ 


Da p D b 

V ' bb ~ V ' 


P h — 


q-y \y7ra yirbj 7rab T yA a yKbJ KabJ • 


- Then 




Ptt-x — 


9-nafTxaPaa + ^g-Kag-Khf-na firbPab 


Tkk = 


Tttk '■ K — > K, 


TttK = 


Tk-k 



(11) 



— g-iragKafira fKaPaa + {g-rragKb fjraj Kb 

^-9Ka9nb fxaf-KbjPab + QitbQKbfitb fxbPbb, 

(12) 

is a set of unitary elastic amplitudes. 

- The inelasticity i] is related to the amplitudes or S- 

matrix as: 



ve 2i5 P = Spp = 1 + 2 ippTpp , P = tt, K 



\/T - V 



2jS 7rK _ 



— — i S^k — Z^/PttPkTkk , 



where the sum of pion and kaon phase shifts is: 
Kk = K + $k ■ 



(13) 



(14) 



In the following, we shall work in the minimal case 
with one shape function: 



fna(s) = fnb(s) = fKa(s) = f Kb{s) , 

where: 

0£> = V~Dtt = &DK ; SA = SAtt = S AK 



(15) 



(16) 



3. Phenomenology of elastic tttt —> tttt scattering 
Data input 

The only data input used in this process is the pion 
phase shift 8^ well measured experimentally. We shall 
use the new precise data from NA48/2 on KeA = K — > 
■KiteVe for the 7T7T-phase shift below 390 MeV [49] and 
use from 400 to 900 MeV the CERN-Munich [50] and 
Hyams et al. |51] 7T7r-phase from tttt — > 7T7t which agree 
each others above 400 MeV. These data are shown in 

Fig. m 

"bare" resonance = A</> 4 model 
Let's first fit the elastic tttt data by using a A</> 4 model 
without any "bare" resonance. In this old version of the 
model [S5J , one can introduce the shape function fp [57] : 



Tpp — 



A/p(a) 



1-A/p(s 
where o~r>\ = &d-k and: 

Us) = -\h 2 (s) - h 2 (0) 

TT L 

with : 

hp{s) = fp(s)L s i(s) 



fp(s) = 



SAP 



(s + o- m )(s + a D2 ) 



E 



0~Ni 
+ 0~Di 



L s i( — 0-Di 



,(17) 



(18) 



(19) 
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Figure 1. Central values of the best fits of the irir phase shift 
below v / s=0.75 GeV for elastic irir —> irir scattering: blue (dotted) 
line for "bare" resonance with x^ jn reached at y/s= 0.7 GeV; 
red (continuous) line for 1 "bare" resonance and X^in reacne d at 
t/s= 0.7 GeV; green (dashed) line for 2 "bare" resonances and 
Xmin reached at v / i= 0.75 GeV. 

&ni,o'N2 in Eq. (H"7|) are the residues of fp(s) at 
<jd\ = o\Dtt, <JD2- In fitting the "bare" parameters, 
we look for a minimum of x 2 = Xmin by varying the 
range of the interval [4m \ , s] inside which we per- 
form the fit. Here, this is obtained for y/s=0.7 GeV 
where: xL„/ndf=12.04/14=0.86. The fitted values of 
the "bare" parameters and the resulting values of the 
physical pole parameters are given in TablefTpl. The 
quoted errors of the "bare" parameters come from the fit 
program MINUIT. The errors induced by each of these 
"bare" parameters on the physical poles can be added 
(as currently done) linearly or quadratically H. These 
results indicate that, though not accurate, this original 
version of the model gives a reasonnable value of the 
physical parameters. 
1 "bare" resonance 

This analysis has been done in [47] using the CGL 
parametrization based on Roy equations with con- 
straints from chiral symmetry |58) . In the following, we 
shall use instead the new precise data from NA48/2 on 
K e4 for the 7T7r-phase below 390 MeV and use from 
400 to 900 MeV the CERN-Munich [50. and Hyams et 
al. [51] 7T7r-phase from tttt — > tttt which agree each oth- 
ers above 400 MeV. We extract the "bare" parameters 
from these data: 

In the first step, we leave all "bare" parameters free 
and find a minimum % 2 : x^ in /ndf=9. 43/17=0. 55 for 
\fs = 0.75 GeV. The fitted value of the Adler zero is: 



SA „ = 0.0394(92) GeV 2 , (20) 

5 Fixing the Adler value at the one in Eq. 121 1 does not bring any 
improvements here. 

"Notice that if we have added linearly the errors induced by each 
"bare" parameters by taking into account their signs, we would 
have obtained about 2 times more accurate predictions. For a 
more conservative error, we shall take here and in the following 
the quadratic sum of these errors. 



Table 1 

Values in GcV d (d — 1,2) of the bare parameters of the K-matrix 
model for different "bare" resonances input for tvtt — > tv-k elastic scat- 
tering. The fit has been performed until y/s ~ 0.7 GeV (0 res. and 1 
res.) and 0.75 GeV (2 res.), where the \ /ndf is minimal. The cor- 
related errors of the "bare" parameters come from the fit of the data 
using MINUIT. The ones of the physical poles arc the quadratic sum 
of the errors induced by each "bare" parameters (a linear sum would 
lead to 2-3 times more accurate values due to cancellations of some of 
the errors in this case.). An average is given in the last column. 



Output res. 



Ires. 



2 res. 



Average 



SA 

O'Dtt 

&D2 

S Ra 

A 

a 

9irb 

xL„ 

ndf 



0.009(6) 
6.2(3.2) 
7.6 ±4.5 

108(34) 



fixed 
1.41(7) 

1.94(9) 

2.54(8) 



12.04 
14 



=0.86 = 0.78 



M a 468(181) 456(19) 
T CT /2 261(211) 265(18) 
\g an+n -\ 2.58(1.31) 2.72(16) 



fixed 
1.78(10) 

26.97(1.54) 

10.42(30) 

0.61(31) 

-0.39(8) 

16 U -' J 

448(18) 
260(19) 
2.58(14) 



452(13) 
259(16) 
2.64(10) 



Table 2 

Mass and 1/2 width in MeV of the a meson in the complex plane. 



Processes 



M a - iT a /2 



Refs. 



This work 
KeA 
KeA 
Average 



7T7T 
7T7T 



7T7T 452(13) -i 259(16) 

mr/KK 448(43) - i 266(43) 
452{12) - i 260{15) 



Others 

7T7T — > 7T7T0Roy( 
7T7T — > 7T7T / KKQ 
J/ljj — > LJ7T7T 
D + — > 7T + 7T~7r + 



)ChPT 441+^ 6 - i 272^5 
Roy 461 ± 15 - i (255 ± 16) 
541 ±39 — i (222 ± 42) 
478±29-i (162 ±46) 



59 
60. 

EH 



which is relatively bad compared with the theoretical 
expectation^ : 



sa-k = -f = °- 0094 Gey2 



(21) 



- Then in the second step, we fix the Adler zero at 
the value in Eq. [21] and deduce the results in Tabled] 
The fit is shown in Fig. [TJ These "bare" parameters 
lead to the physical poles in Table [T] which we con- 
sider as improvements of the previous results in [4 7) . 
This result is comparable in size and errors with the 
precise determinations from recent analyses of the anal- 
ogous 7T7T — > 7T7T / KK scatterings data using different ap- 



7 Here and in the following, we shall not explictly include isospin 
breaking effects which we expect to give small corrections. 
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proaches (Roy equations chiral symmetry constraints 
[55] , Roy equations © control of the high-energy be- 
haviour of the amplitude [55]) (Table HJlI, which have 
been obtained before the last KeA NA48/2 precise data 

ESI- 

2 "bare" resonances 

We repeat the previous analysis by working instead 
with 2 "bare" resonances. We fix the Adler zero at 
the value in Eq. [21] and fit the other "bare" param- 
eters. We obtain the results quoted in Table [1] for a 
xL„/ndf=12. 71/16=0. 794 at y/s = 0.75 GeV. The fit 
is shown in Fig. [T] 

Comments and final _results from tttt — > tttt 
From previous studies, we conclude that: 
- The results from different forms of the model in Table 
[T]are very stable. The final results from elastic tttt — > tttt 
are the average of the ones from 0, 1 and 2 "bare" input 
resonances quoted in this Table [TJ which are: 



M CT [MeV] 

IffcTTT+TT- | 



452(13) -i 259(16) 
2.64(10) GeV . 



(22) 



- The results, from the "bare" resonance or \<fi 4 model 
show that the existence of the a pole is not an artifact 
of the "bare" resonance entering in the parametrization 
of the tttt amplitude Tpp. 

Noting that the concavity of the fit curve in Fig. [1] 
around the p-meson mass region has raised some doubts 
on the data of the phase shift 6^ [52], we have redone 
the fit by assuming that the data increases linearly from 
the KeA one. Using 1 or 2 resonances, we still find, in 
this extreme case, a a pole : 



M a [Me V] ~ 413 -i 300 



(23) 



where a similar value has been obtained earlier |62| . 
This result may indicate that the existence and the dy- 
namics of the a is mainly due to the low-energy be- 
haviour of the tttt phase shift 8^ data, which are accu- 
rately determined from KeA by NA48/2 [49] . 

4. Phenomenology of inelastic tttt —> tttt/KK 
2 "bare" resonances 2 channels parameters 

In so doing, we take in Tabled three representatives 
sets of tttt — > tttt/KK data in the existing literature: 



Table 3 

Different data used for each different sets for determining the "bare" 
parameters in Table [4] (S^is the tttt phase shift, 77 is the inelasticity 
and 1S the sum of the it and K phases. 



Input 


Set 1 


Set 2 


Set 3 


8 V 


[49 50 51 


[49150151] 


[49150)51] 


V 


m 


M 


[52] 


8ttK 


53 


52 


m 



NA48/2 
CERN-Munich 
Set 1 
Set 2 
-Set 3 





s"(GeV) 




1.10 115 

s 1,! (GeV) 



Cohen 
Etkin 
Martin 
Set1 
Set 2 
-Set 3 



3 See also Table [2] for some other determinations. 



Figure 2. a) Fit of the 7T7T phase 5^ versus ^/s: Set 1 (blue: 
dotted); Set 2 (green: dashed); Set 3 (red: continuous); b) Fit of 
the inelasticity rj; c) Fit of the sum of it and K phase S^pc ■ 



- The choice for 5^ which we have used in the analysis 
of elastic tttt — > tttt scattering is unique and comes from 
the new KeA data of NA48/2 [49] below 390 MeV and 
from tttt — > tttt/KK data above 400 MeV measured by 
CERN-Munich [SD] and Hyams et al. [5T] [see Figs. □ 
and [2] a)]. 

- For the inelasticity 77, different data exhibits a min- 
imum rjmin just above the KK threshold. CERN- 
Munich and Hyams et al. give the smallest value: 
Vmin ~ 0.4, while Cohen et al. [52] provide the largest 
one: r] min w 0.7 [see Fig. (2]b)]. 

- For the sum of tt and K phase we use the one 
from Cohen et al and from Etkin-Martin [53] , which rep- 
resent the two extreme cases [see Fig. [2]c)]. 

With these choices, we expect to span all possible re- 
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Set 1 
Set 2 
Set 3 



1.05 1.10 1-15 1.20 1.25 1.30 1.35 1.40 

s'°(GeV) 



Figure 3. x 2 / n df versus y/s in GeV : Set 1 (blue: square); Set 2 
(green: circle); Set 3 (red: triangle). 



Table 5 

a and /o(980) meson parameters from tttt —> tttt/KK scatterings using 
the bare parameters in Table HI: the mass and width arc in McV, while 
the couplings are in GeV. The errors are the quadratic sum of the ones 
induced by the "bare" parameters in Tablc[^(a linear sum would lead 
to 2-3 times more accurate values due to cancellations of some of the 
errors in this case). An average is given in the last column. 





Ciot 1 

oei i 


Clot 9 

oet z 


Clot ^ 

oet o 




M a 


435(74) 


452(72) 


457(76) 


448(43) 


r CT /2 


271(92) 


266(65) 


263(72) 


266(43) 


|5(77T+7r~ | 


2.72(78) 


2.74(61) 


2.73(61) 


2.73(38) 


\9aK+K~ 


1.83(86) 


0.80(55) 


0.99(68) 


1.06(38) 


M f 


989(80) 


982(47) 


976(60) 


981(34) 


17/2 


20(32) 


18(16) 


18(18) 


18(11) 




1.33(72) 


1.22(60) 


1.12(31) 


1.17(26) 


\9fK+K~ 


3.21(1.70) 


2.98(70) 


3.06(1.07) 


3.03(55) 



gions of the space of parameters, and then to extract re- 
sults which do not only come from a single experiment. 
We have not used the data of Kaminski et al. {54] due to 
the large errors, which, however, agree within the errors 
with the other data sets used here. 



Table 4 

Values in GcV d (d — 1,2) of the bare parameters of the K-matrix 
model for 2 channels © 2 bare resonances from Ke4 ® tttt — y tttt/KK 
scatterings. The fit has been performed until y/s ~ 1.225 ~ 1.250 
GeV, where the /ndi is minimal (see Fig. [3j. The correlated errors 
come from the hitting procedure using the program MINUIT. 



Output 


Set 1 


Set 2 


Set 3 


SA 


0.016 ±0.004 


0.013 ±0.006 


0.010 ±0.006 




0.740 ± 0.097 


0.909 ±0.201 


1.116± 0.262 


S Ra 


4.112 ±0.499 


2.230 ±0.271 


2.447 ±0.298 


9na 


-0.557 T 0.177 


0.846 ±0.391 


0.997 ±0.516 


gKa 


3.191 ±0.499 


1.458 ±0.262 


1.684 ±0.363 


S R„ 


1.291 ±0.062 


1.187 ±0.094 


1.354 ±0.149 




-1.562 t 0.117 


-1.527t0.134 


-1.756 T 0.183 


gxb 


0.748 ±0.062 


0.999 ± 0.149 


1.159 ±0.261 


v 2 . 
V-m i n 

ndf 


^=0.914 


^ = 0.759 


^ = 0.763 



In the following, we shall use: 

1 

2 



m K 



(m K + + m K o) = 495.65 MeV . 



(24) 



Letting all "bare" parameters free, we study in Fig. [31 
using the fitting program MINUIT, the variation of 
X 2 /ndf versus t/s until 1.4 GeV where the data are avail- 
able. In the fitting procedure, we have chosen the same 
initial conditions for the 3 sets, where a good conver- 
gence with a good x 2 / n df of the solutions has been ob- 
tained for Set 2 and Set 3. A minimum value Xmi„/ndf 
is reached for ^fs ~ (1.225 — 1.250) GeV at which we 
extract the optimal outputs given in Table [4] At each 



corresponding value of Xmm / n< if, the fits for different 
sets of data are shown in Fig.[2J All three sets give good 
values of x^ in /ndf less than one. 
Poles from 2 "bare" resonances 2 channels 
We use the results of the "bare" parameters in Table@] 
obtained at Xmin/ n ^ f° r deducing the ones of the com- 
plex poles in Table[5] The errors on the physical poles 
are induced by the ones of the "bare" parameters in 
Table|4] and have been added quadratically. The itera- 
tion of solutions from Set 1 has only a local minimum 
in x 2 such that, in order to be more conservative, we 
have multiplied by a factor 2 the related uncertainties 
of the results coming from the fit. The last column 
gives the mean value from the three different determi- 
nations. We have taken (as is usual in the literature) 
the weighted average, where the corresponding error is 
more weighted by the most accurate predictions^. One 
can see in Table[S] that the results from different sets 
of data are unexpectedly stable for both a and /o(980) 
parameters, which increase our confidence on their in- 
dependence on the input data sets. 
- For the a, we obtain the average of the complex pole 
mass and width given in Table [5J 



M a [MeV] = 448(43) - i 266(43) 



(25) 



This result is in perfect agreement with the mean value 
from elastic 7T7t — > tttt scattering in Eq. [22] and compara- 
ble in size and errors with the ones in TableO Averaging 
the two predictions in Eqs. [22] and [25] we deduce our 
final value: 



M a [MeV] = 452(12) - i 260(15) . 
Averaging the result in Table [5] and Eq. 
Ifc+rl = 2.65(10) GeV, 

= 0.37(6) , 



\9aK+K- I 



\g a 



(26) 
we deduce: 

(27) 



9 Alternatively, we can take, as a final value, the most accurate 
prediction, which leads about the same result. 
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which improves and confirms our previous rough find- 
ings in [?5] and which is comparable with some other 
determinations in Table[5]from 48 0: the sizeable cou- 
pling of the a to KK disfavours the usual tttt molecule 
and four- quark assignement of the a, where this coupling 
is expected to be negligible. 
- For the /o(980), we obtain the mean value: 

M/[MeV] = 981(34) - i 18(11) , (28) 

which is comparable with the PDG range of values [64] : 

M/pVIeV] = 980(10) - i(20 ~ 50) . (29) 

From Table [SJ we also find: 



\g fw+w -\ = 1.12(31) GeV 

_ \9fK+K- 



WfTT + TT- | 



2.59(1.34) , 



(30) 



in agreement with the determinations in the existing 
literature (see Table[6]). The large value of this ratio 
of coupling and the relative narrowness of the /o(980) 
width (compared to e.g. the p-meson) does not favour 
the pure (uu + dd) content of the Jo (980) where r^K 
is expected to be about 1/2 and the width of about 120 
MeV [6 91. This feature has been used as an indication 
of the four-quark nature of the Jo (980) (see e.g. 163]) 
or alternatively of its large gluonium component via a 
maximal mixing with a qq state (see e.g. f5!6] /). 



Table 6 

Modulus of the 7r^7r~ and K~ complex couplings in GeV of the 
(7 and of /o(9S0) from S- and K-matrix models for tttt — > tttz/KK 
scatterings compared with the ones from <fi and J/ip decays. rsnK = 
\9sk+k- \/\9 s *+*- l : S = CT > /■ 



Processes 



Igq-TT+TT- I r*KK Iff/TT + TT-I ModdS 



This work 

7T7T -> rnr/KK 2.65(10) 0.37(6) 1.17(26)2.6(1.3) [47155] 
Others 

7T7T irir/KK 2.03(3) 0.65(18) 0.97(6) 1.7(2) [54] 

2.5 0.62 1.55 1.20 [65] 

0-> cr// (980) 7 - 0.67 [66] 

rp -> (f> nn/KK 2.35 1.80 [57] 



Average 



2.4 



0.6 



1.5 



Model with 1 "bare" resonance © 2 channels 

We have further studied the influence of some other con- 
figurations of the model on the fit of the a by analyzing 
the minimal case: 1 "bare" resonance © 2 channels and 
using for instance Set 3 of data. Letting all "bare" pa- 
rameters free, the Xmin/ n ^ — 10 is very bad which is 
obtained by doing the fit from Im^ to s/s = 1.2 GeV. 



10 Similar values of r a ^K are also found from some fits in |63| but 
the results obtained there are unstable. 



As (intuitively) expected, the previous results for the a 
parameters are approximately reproduced: 



M a [MeV] ps 377 -i 195 , 
and 

l<?/7r+7r-| ~ 2.13 GeV, r^A-wO.42 



(31) 



(32) 



while the fo mass is pushed far away from the KK 
threshold: 



M/[GeV] «3.8 + i 1.7 



(33) 



Due to the bad quality of x 2 / n df, the result from this 
version of the model will not be retained. 

5. On-shell mass, width and couplings of the a 

Due to the large width of the c, a direct comparison of 
the previous results with the ones obtained from QSSR 
or some other theoretical predictions in the real axis is 
questionable. For better comparing the results obtained 
in the complex plane with the theoretical predictions 
obtained in the real axis, it is more appropriate to in- 
troduce like in [ITJ the on-shell meson [55] masses and 
hadronic widths, where the amplitude is purely imagi- 
nary at the phase 90°: 



Rc2?((M° s ) 2 ) = 0=^M° S ~ 0.9 GeV 



(34) 



In the same way as for the mass, one can also define an 
"on-shell width" [47] from Eqs. © and (0 evaluated at 



a = (Af° s ) 2 : 



^ " Re^ 



0.7 GeV 



(35) 



which are comparable with the Breit-Wigner mass and 
width [51169170] : 

Mbw — ^ bw — 1 GeV . (36) 
These values lead to the on-shell coupling: 
|<C +7r _|~6GeV. (37) 

6. Comparison with QSSR © LET predictions 

- One on hand, the corresponding on-shell (or Breit- 
Wigner) mass and coupling of the a are comparable in 
size with the predictions from combined QSSR © LET 
analysis [6 30 32. for a glueball with a large OZI viola- 
tion for its coupling to tttt and if .Km : 

M a ~ 1 GeV , \ 9(Tn+n - 1 ~ \g aK+K - 1 ^ 5 GeV , (38) 

implying: 



a— *-7T + 7T - 



VcmM n 



1- 



4m 2 



0.7 GeV 



(39) 



The existence of the a is necessary for a consistency 
between the subtracted and unsubtracted QSSR [32] . 



11 The expectation of a glueball chiral coupling to pair of Gold- 
stone bosons |24| could not hold in this non-perturbative regime. 
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where the gluonium two-point correlator subtraction 
constant [35] : 

m ^-^(a s G 2 ) (40) 

16 7T 

plays a crucial role = — 11/2 + n/3 for n flavours), 
and where the value of the gluon condensate is: 
(a s G 2 ) = (6.8 ± 1.3) x lO" 2 GeV 4 [71l72l73l74] , 

- On the other hand, QSSR predicts for & S2 = 
l/^/2(uu + dd) 1=0 scalar meson [619j: 

M S2 ~ 1 GeV , T S2 ~ 0.12 GeV , (41) 

and 

|9s 2 7r+7r- 1 — 2.5 GeV , \f*±*-\ = 1 (42) 
These results indicate that: 

-The 5*2 is narrower and much higher in mass than the 
complex a pole often identified with a qq state in the 
current literature. 

- The /o(980) cannot be a pure (uu + dd) state due to 
the large ratio of its KK over its nn couplings rf„K 
[Eq.[30] and to its relative (compared to the p-meson) 
small 7T7T width. It cannot be also a pure ss or KK 
molecule due to its non-negligible width into 7T7t. 

- A large gluonium component eventually mixed with a 
qq state in the a and / (980) |2l5l6l22l32l47l48j can be 
advocated for evading these previous difficulties. 

- Further phcnomcnological searches for gluonium have 
been proposed in the literature in <f>, J / ip and T radiative 
decays [BI9I32] . D, B semi-leptonic [75 9 and hadronic 
[22] decays. 



r CT7r x ~ 0.37(6) [Eq. [27] is a strong indication against a 
pure 7f7r molecule and four-quark substructure of the 
<7, whilst its large width cannot be explained (using 
the QSSR results in a previous section) from a simple 
(uu + dd) assignement. 

- The large value: r f „ K ~ 2.59(1.34) [Eq. [30] of the 
ratio of the /o (980) couplings to KK over the one to nir 
and of the /o(980) relative narrow width compared e.g. 
with the one of the p-meson does not favour the pure 
(uu + dd) assignement of the /o(980). In this scheme 
one would predict a ratio of coupling of about 1/2 and 
a width of about 120 MeV [Eq. Hi] . 
The four-quark scenario can explain the large KK cou- 
pling of the fo (980) but it fails to explain the large cou- 
pling of the a to KK. 

The simple qq scheme cannot explain the large KK cou- 
pling and the narrowness of the /o(980) as well as the 
broad width of the a. 

A large gluonium component eventually mixed with 
a qq state of the a and /o(980) can be advocated 
[2|5l6l22|32l47|48j for evading these above-mentioned 
difficulties. 
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7. Summary and conclusions 

We have used new KeA = K — ¥ ~K-Kev e on rnr phase shift 
below 390 MeV © different tttt — »• ttit/KK scatterings 
data above 400 MeV, for extracting the a = Jo (600) and 
/o(980) masses, widths and hadronic couplings, within 
an improved "analytic if-matrix model" . 
Using a A0 4 version of the model, we have noticed from 
different analysis that the existence of the a in the com- 
plex plane and having a mass of about 452 MeV is not 
an artifact of "bare" resonances used in the analytic K- 
matrix model. 

We have also seen that our predictions are very stable 
vesus the different forms (number of "bare" resonances) 
of the models. Our results are summarized in Table[5] 
The masses and widths [Table [2] and Eqs. [26] of the 
a are comparable in size and errors with the most ac- 
curate determinations in the existing literature [58 59 
(see also Table [2]), while the ones of the /o(980) in Eq. 
[28] are comparable with the PDG values [64]. The small 
uncertainties in our determinations can be mainly due 
to the new accurate data on Ke4 from NA48/2. 
The values of the couplings confirm and improve our pre- 
vious results in [48] and are comparable with the ones 
from some other processes given in Table [6] 
- The (unexpected) sizeable coupling of the a to KK: 
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